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TRANSPORT  IN  THERMALLY  GROWN  SILICA  ON  SILICON 


INTRODUCTION 

Computer  modeling  of  the  oxidation  stop  in  silicon  integrated  circuit 
(IC)  device  processing  requires  an  improved  understanding  of  the  nature  of 
oxidant  transport  and  the  oxidation  mechanism,  on  a  microscopic  scale,  so 
that  an  accurate  mathematical  description  ran  he  formulated.  The  oxidation 
of  silicon  has  developed  along  three  almost  separate  paths.  The  first  Is  thg 
very  high  temperature  oxidation  from  about  llOD'V  to  above  the  melting  point, 
1420°C,  where  the  oxidation  may  involve  volatile  SiO  gaseous  species  and 
transport  through  a  boundary  layer  may  he  rate  determining.'-^  The  second 
is  the  high  temperature  oxidation  from  about  900  to  1  TO  )°C  in  which  a  thick 
layer  of  silica  separates  the  reactants  and  a  Jccrease  in  rate  with  time  is 
observed.  The  studies  in  this  regime  have  generally  been  directed  to  protec¬ 
tion  of  ceramics  and  metals. Frequently  the  material  studied  has  been  im¬ 
pure  silicon  or  silicon  carbide  or  silicon  nitride.  The  kinetics  have  gener¬ 
ally  been  described  as  parabolic, 

2 

x  =  k^.t 

where  x  denotes  film  thickness,  t  the  time, and  kj  is  the  Tamman  parabolic 
rate  constant  expressed  in  units  of  cm^/sec.  The  most  extensive  studies  have 
been  carried  out  at  lower  temperatures,  about  500  to  1000°C,  on  relatively 
pure,  well-oriented  single  crystals  for  the  electronics  industry.  In  this 
regime,  the  mixed  linear-parabolic  kinetic  model  first  described  by 
U.R.  Evans^t  and  applied  especially  to  silicon  by  Deal  and  Grove®-?-  is  supposed 
to  hold.  The  advantage  of  the  single-crystal,  lower-temperature  studies  is 
that  the  materials  are  well-characterized  and  controlled.  The  advantage  of 
the  higher  temperature  studies  is  that  the  role  of  minor  impurities  should 
be  of  less  influence  because  of  the  thermal  disorder  in  the  grown  silica 
layer.  The  very  high  temperature  work  will  not  be  considered  further  here 
inasmuch  as  transport  through  a  solid  silica  layer  is  the  aim  of  this  report. 

KINETIC  EQUATIONS 

The  equation  for  the  mixed  linear-parabolic  model  due  to  Deal  and 
Grove®^  is  in  integrated  form, 

x^  +  Ax  =  B(t  +  t) 

where 

A  =  2Deff[_l/k  +  I/hJ 

B  -  2DeffC*/Nl 

T  =  (x^  +  Ax^) /B 

x  =  thickness  of  Si02  at  any  time,  t 

=  thickness  of  oxide  on  Si  initially  present  prior  to 
onset  of  kinetic  study,  x  =  at  t  =  0 

D^jj.  =  effective  diffusion  coefficient  for  species 
migrating  through  Si02 


t  See  list  of  references  on  page  12. 
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k  =  r.itt-  consign* 


h  =  rate  constant 

C*  =  concentration  of  oxvgen  at  the  0,-Sii),  interface 

Nj  =  number  of  oxidant  molecules  incorporated  into  a  unit 
volume  of  oxide  laver. 


In  the  limit  of  thick  films  (x 
the  simple  parabolic  rate  expression 


)  this  equation  should  approach 


=  Bt  =  kTt 


4  5  7-19+ 

The  oxidation  kinetics  have  been  studied  by  many  investigators.  *  ’ 

Most  of  the  data  have  been  obtained  at  one  or  near  one  atmosphere  oxygen 
pressure.  When  parabolic  kinetics  are  observed,  generally  diffusion  of  the 
species  across  the  reaction  product  layer  is  rate  determining.  Under  such 
conditions  and  when  there  is  virtual  equilibrium  at  the  metal-oxide  interface 
and  at  the  oxide-gas  interface,  the  classical  equation  for  oxidation  of  metals 
by  C.  Wagner20,21,  may  appiy.  Using  the  Tamman  rate  expression 


x  =  kTt, 

it  can  be  shown  that. 


where 


(a  )  a 
ion  elec 


k  =  rational  rate  constant  in  eq/cm  sec 
r  n 

v  =  equivalent  volume  (eq  wt: density) 

Nq  =  Avogadro's  number 
F  =  Faraday  constant 


o  =  total  conductivity  =  a.  +  o  , 

ion  elec 

Uq  =  chemical  potential  of  nonmetal  (in  this  case  oxygen)  at 
the  metal-oxide  interface  (Si-SiO^  interface) 

Uq  *  chemical  potential  of  oxygen  at  the  SiC^-gas  interface 

Zg  =  valence  of  oxygen  in  SiC^ 


+  See  list  of  references  on  page  12. 
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cion  »  ionic  conductivity  =  in  Si00 

a  .  =  electronic  conductivitv  in  SiO,. 

elec  •  2 

In  order  to  test  the  applicability  of  the  C.  Wagner  model,  information 

is  needed  concerning  the  values  of  c.  and  J  ,  as  a  function  of  oxvgen 

,  ion  elec  • ^ 

potential,  jjq  =  +  RT  ln  P0^  • 

Three  types  of  data  are  needed.  These  are 

1.  isothermal,  isobaric  oxidation  kinetics  to  test  whether  the 
parabolic,  mixed  linear  parabolic  or  other  rate  equation  holds. 

2.  isothermal  oxidation  kinetics  as  a  function  of  oxygen  partial 
pressure.  These  data  can  be  usgd  to  test  an  oxidation  mechanism 
according  to  Fueki  and  Wagner?^ 

3.  isothermal  ionic  and  electronic  conductivity  of  Si02  as  a  function 
of  chemical  potential  of  oxygen. 

ELECTRICAL  CONDUCTIVITY  AND  IDENTIFICATION  OF  MIGRATING  SPECIES 

Additionally^,  the  diffusing  species  should  be  identified.  Marker 
studies^* 23-26 .  have  indicated  that  the  migrating  species  through  the 
silica  layer  is  most  likely  oxygen  and  not  silicon.  The  charge  on  the 
oxygen  and  whether  it  is  a  single  atom  or  a  molecule  is  a  matter  of 
controversy. 

61“ 

Deal  and  Grove  '  suggested  that  the  migrating  species  was  the  peroxide 
ion,  O^.  Others  have  suggested  that  the  oxygen  molecule  is  uncharged  or 
that  an  interstitial  oxygen  ion,  0 is  the  dominant  species.  These  cases 
will  be  discussed  further  below.  Jorgensen^»27-29+  carrieij  0ut  a  series  of 
classical  papers  using  an  applied  electric  field  to  accelerate  and  to  retard 
(depending  on  the  polarity)  oxidation  of  silicon  and  of  other  metals.  The 
fact  that  an  applied  electric  field  could  accelerate  or  retard  the  oxidation 
rate  was  accepted  by  many  authors  as  evidence  that  the  oxygen  species  pos¬ 
sessed  a  negative  charge  or  charges.  However,  Raleigh^O  pointed  out  that 
this  was  not  necessarily  the  case.  An  applied  electric  field  could  result 
in  electrolysis  of  the  oxide  film  and  the  existence  of  a  charged  species 
could  not  be  unambiguously  determined  from  such  experiments  alone.  Mills 
and  Kroger-*1'  carried  out  experiments  on  silica  to  obtain  the  total,  ionic 
and  electronic  conductivity.  They  used  the  following  configurations: 


Cell  I 

1  S1 

Si02 

Si 

Cell  II  1 

Si 

[Si02 

!  Gei 

Si 

-x  x 

Cell  III  J 

Si 

Si02 

-d 

O 

ho 

i 

t  See  list  of  references  on  page  12. 


Cell  I  is  a  symmetrical  cell  in  which  the  chemical  potential  of  silicon 
has  its  lowest  value  (ag^  =  1).  Current-voltage  characteristics  were  used 
to  obtain  the  total  conductivity,  c(ionic  plus  electronic;  c’  =  ~'ion  +  'elec) 
using  this  cell.  The  results  from  this  particular  cell  have  been  criticized 
on  the  basis  that  the  dielectric  breakdown  of  silicon  was  exceeded.  An 
approximate  check  on  this  can.be  made  as  follows.  The  dielectric  strength 
of  Si02  is  "c  5x10^  volts/cm^  presumably  at  25°C.  The  thickness  of  Mills 
and  Kroger's  films  was  3000  %  and  the  applied  voltage  was  0.1  to  10  volts. 
Hence  the  field  they  applied  was 


V  VO  1  t  q  —  Q  A 

—  in  -  -  (0.1  to  10)  •  3000x10  =  3x10  to  3x10  volts/cm. 

x  cm 

The  temperature  rr-nge  Mills  and  Kroger  used  was  150°C  and  above.  This  simple 
calculation  indicates  their  measurements  at  the  elevated  temperatures  ap¬ 
proached  the  dielectric  breakdown  strength  and  hence  further  tests  of  this 
cell  using  different  oxide  thicknesses  would  be  highly  desirable. 

Cell  II  is  an  asymmetrical  cell  for  measuring  the  open  circuit  potential, 
E,  (the  en.f  at  infinite  impedance)  across  a  layer  of  Si02-  This  voltage  is 
related  to  the  oxygen  pressure  by 


P02’  P02  =  oxygen  partial  pressure  at  the  Si~Si02  and  at  the  SiC^*- 
Gej_xSi^  phase  boundaries,  respectively 

t.  =  transport  number  for  ions 
ion 

n  =  number  of  equivalents  transported  for  a  reaction  invol/ing 
1  gm  mole  C^.  If  0_  is  the  migrating  species  then  n=4 . 

Thus  by  varying  the  p02  (changing  the  value  of  x  in  the  alloy  Ge^^Si^) 

the  isothermal  ionic  transport  number  as  a  function  of  oxygen  pressure  or 
Si/0  ratio  can  be  obtained. 


Cell  III  is  an  analogous  asymmetric  cell  to  Cell  II  and  the  same  prin¬ 
ciples  apply. 

Therefore  if  the  total  conductivity  can  be  obtained  from  Cell  I  and  the 
ionic  transport  number  from  Cell  TI  or  Cell  III,  then  the  ionic  conductivity, 
O^on,  may  be  obtained  as 


and  the  information  used  to  test  the  applicability  of  C.  Wagner's  equation 
plus  an  indication  of  the  state  of  charge  on  the  migrating  oxygen  (the  value 
of  n  in  the  equation  above). 

t  See  list  of  references  on  page  12. 
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Three  sets  of  data  are  available  'or  this  purpose:  the  data  of  Mills 
and  Kroger^!-1-  cited  above,  the  analysis  of  Jorgensen's  data  by  Raleigh^®’1’ 
and  a  report  by  Tripp  et  al.^+  The  problem  is  that  if  one  extrapolates  these 
data  to  a  common  temperature  of  850°C,  the  data  of  Mills  and  Kroger  and  of 
Tripp  et  al  differ  by  one  order  of  magnitude  and  that  of  Jorgensen  by  four 
or  five  orders  of  magnitude.  Furthermore,  none  of  the  measured  conductivity 
data  allow  even  an  order  of  magnitude  calculation  of  the  rate  of  oxidation. 

In  addition,  the  state  of  charge  on  the  oxygen  atom  or  molecule  is  not 
clear . ^ H" 

ANALYSIS  OF  OXIDATION  KINETICS 

A  reevaluation  of  the  basic  oxidation  kinetics  is  needed.  Deal  and 
Grove, ^  Ligenza.^b  Hess  and  Deal, Hi"  Attala,^"  Laverty  and  Ryan,^4' 

Tripp  et  al,^*  Flint, ®J"  and  Nakayama  and  Collins-H+  all  report  parabolic 
or  mixed  linear  parabolic  kinetics  to  best  represent  the  oxidation  data. 


The  oxygen  pressure  dependence  of  the  oxidation  kinetics  may  also  be 
used  to  test  the  mechanism  and  the  charge  on  the  diffusing  species.  Table  1 
lists  available  information  concerning  the  oxygen  pressure  dependence  of 
the  rate.  For  the  most  part  the  rate  is  proportional  to  Pq  where 

0.75  <_  m  1.  Recently,  Zeto^i"  completed  a  study  of  the  dry  high  pressure 

oxidation  of  silicon  (up  to  500  atmospheres)  and  showed  that 

1/0.7  1.4286 

x  =  x  -  kzt 

for  isothermal,  isobaric  oxidation.  Furthermore 
k  «  P  °*7 

Z  o2 

The  lattice  defects  in  Si02  can  be  formally  represented  as  point  defects 
in  an  ionic  compound.  These  are  listed  in  Table  2.  The  nonstoichlometric 
defects  corresponding  to  simple  point  defects  are  listed  in  Table  3. 

1 8t 

Irene  very  graciously  supplied  Zeto  a  copy  of  his  data,  film  thick¬ 
ness  versus  time,  for  three  temperatures,  780,  893  and  980°C,  for  oxidation 
of  (100)  silicon  wafers  in  dry  (<  I  ppm  H20(g)  oxygen  (Pq2  =  1  atm). 

181* 

Irene  and  van  der  Meulen  reported  that  the  following  equatign  well 
represented  the  data  beyond  an  initial  thickness,  xq(50  Xq  £  380  A) 
depending  on  the  temperature: 

(t-tQ)  =  A' (x-Xq)  +  B'(x2-Xq) 

where  A’  =  reciprocal  of  the  linear  rate  constant 

B'  =  reciprocal  of  the  parabolic  rate  constant 


Xq, tg  =  film  thickness  and  time  prior  to  the  onset  of  the 
mixed  linear  parabolic  model 


t  See  list  of  references  on  page  12. 
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Dependence  of  oxidation  rate  on  oxygen  partial  pressure 


Formal  representation  of  possible  lattice  disorder  In 


TABLE  3 


Formal  representation  of  nonstoii  hiometric 
point  defects  in  SiO, 


Compensating 
Interstitial  electronic 


species 


species 


Electroneutrality 

condition 


Oxygen  pressure 
dependence 


1 .  uncharged  O, 

none 

— 

po, 

2. 

u ; 

h" 

[o/  ;  *  [  h°  ] 

po/2 

3. 

or 

2h° 

[or  ]  =  i  [h*  ] 

< 

4. 

o ' 

h° 

[o']  =  [h“  j 

po/4 

5. 

o" 

2h° 

[O"]--  J[h*] 

Pr!/6 

Note:  There  is  the  possibility  that  the  concentrations  of  electronic  defects,  in 
these  examples  electron  holes. are  much  less  than  the  concentration  of 
lattice  defects.  For  example,  in  example  2,  [O/]  ^  constant  »  [h°  ]. 
Hence  the  rate  will  depend  on  ttie  diffusion  of  the  electron  holes  and  [h°  ]  uc  p  ' 
Likewise  for  example  3,  [h°  ]  oc  pn‘  /2  ;  for  example  4,  fh’  j  oc  pfl  ‘  1  and  forUi 

example  5,  [h“]  oc  p'/4  .  2  '  ' 

^2 
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On  the  other  hand,  the  data  were  reanalyzed  at  Fort  Monmouth  and  it  was 
found  that  the  simple  equation 


1/0.7  ,  „ 

x  =  kzc 


in  every  instance  represented  the  data  to  a  greater  degree  et  confidence  on 
a  least  square  analysis  than  did  the  mixed  linear  parabolic  model  of  Deal 
and  Grove6,  or  of  Irene  and  van  der  Meulen.18"  Moreover,  use  of  the  simple 
equation  did  not  involve  the  discarding  of  anv  of  Lhe  data  points.  In  other 
words,  all  data  points  fitted  the  simple  equation  with  greater  ont idence 
limits  than  did  the  mixed  linear-parabolic  which  involved  the  discarding  of 
the  first  data  points  in  the  thin  film  regime.  Zeto^*  had  found  that  all 
of  his  data,  including  thick  film  data  (up  to  a  micron  and  thicker)  fitted 
the  simple  equation  using  all  data  points  with  much  better  confidence  limits 
than  any  of  the  mixed  linear-parabolic  models. 


Unfortunately,  many  of  the  literature  references  do  rot  contain  the  raw 
data  even  in  the  form  of  small  scale  figures  so  a  direct  comparison  with 
other  data  is  not  possible  at  this  time.  Another  vet;,  important  point  is 
the  role  of  water  vapor  on  the  oxidation  kinetics.  It  is  well-known  that  the 
presence  of  water  vapor  lowers  the  activation  energy  and  increases  the  oxida¬ 
tion  rate.  But  it  was  only  recently  determined  that  trace  amounts  of  water 
vapor  exert  a  large  influence  on  the  kinetics.  Irene' 7X  has  shown  very  con¬ 
clusively  that  as  little  as  25  ppm  water  vapor  in  one  atmosphere  of  oxygen 
can  affect  the  oxidation  kinetics. 


As  a  result  of  Irene's  studies  as  well  as  those  of  other  researchers 
industry  has  adopted  stringent  requirements  for  the  presence  of  water  vapor. 
However,  much  of  the  data  published  prior  to  about  1975  must  be  examined 
critically  because  it  was  not  known  that  traces  of  water  influenced  the 
kinetics  to  the  extent  demonstrated  by  Irene.  Tripp  et  al^"'  were  aware  of 
this  problem  and  their  high  temperature  ( 1200-1400°C)  data  follow  the  para¬ 
bolic  rate  law,  x2=kxt,  very  well.  One  of  these  workers,  H.C.  Graham,  pro¬ 
vided  copies  of  that  group's  raw  data.  Their  data  did  not  follow  the 
equation. 


Thus  a  most  pressing  need  is  for  oxidation  kinetics  between  about  700  and 
1200°C  in  dry  (<  1  ppm  water  vapor)  oxygen. 


The  microstructure  of  the  oxide  should  be  examined  as  a  function  of  time 
and  temperature  as  well  as  oxygen  pressure  using  optical  microscopy,  scanning 
electron  microscopy  (SEM),  and  transmission  electron  microscopy  (TEM) .  The 
conditions  under  which  the  two  rate  equations,  x-^k^t  and  x'/0-7=kzt,  hold 
must  be  settled.  Concurrently,  because  the  oxygen  pressure  dependence  is 
not  in  accord  with  the  simple  defect  models  (see  Tables  1  and  3),  and  because 
amorphous  silica  is  a  polymer  or  a  glass  (35-37)  alternative  models  for 
transport  should  be  sought.  For  example,  polymeric  materials  may  exhibit 
defects  due  to  chain  folds.  Likewise  the  presence  of  micropores  or  channels 
must  be  tested. 38+  Silica  is  very  nearly  stoichiometric.  The  excess  oxygen, 
i.e.,  the  oxygen  in  the  silica  layer  at  the  silica-oxygen  interface,  has  been 
reported6''’  to  be  5.2xl0'6/CC  at  1000°C.  This  corresponds  to  about  0.0001  at.Z 


t  See  list  of  references  on  page  12 


excess  oxvgen.  Consequently,  sna 1 1  amounts  or"  impurities  mav  cause  large 
changes  in  the  rates  of  oxidation.  Revesz  and  Evans  demonstrated  this 

effect  dramatically  by  utilizing  different  reaction  tubes,  tvpes  of  heating 
and  sodium  content.  It  is  known  that,  e.g.,  traces  of  iron  catalyze  the 
crystallization  of  amorphous  silica  to  cr is tobal i te . **"  The  activation  ener¬ 
gies  for  oxidation  when  cristobalite  forms  increase  but  the  rates  are  more 
rapid."*  This  indicates  the  preexponential  factor  in  the  Arrhenius  equation 
must  also  be  changing  dramatically.  Again,  a  study  of  the  microstructure 
and  microchemical  analysis  of  the  silica  layer  is  much  needed.  Also,  a  study 
of  the  effect  of  n  and  p  substrates  and  an  analysis  of  the  incorporation  of 
the  dopants  in  the  silica  layer  are  needed.  Here  an  ion  milling  unit  and 
Auger  analysis  would  be  helpful.  A  summary  of  some  of  the  activation  ener¬ 
gies  reported  for  supposedly  amorphous  silica  grown  in  supposedly  dry  oxygen 
is  shown  in  Table  U.  The  energies  fall  into  three  categories.  These  are 
65  kcal  (due  to  formation  of  cristobalite);  '  -10  kcal  (in  dry  oxygen);  and 
30  kcal  or  less  (in  moist  oxvgen). 

SUMMARY 


Based  on  the  recent  data  of  Zeto  and  analyses  of  literature  data,  there 
is  a  regime  of  oxidation  kinetics  which  obeys  the  equation. 


1/0.7 

x 


V 


with  greater  confidence  limits  than  the  simple  parabolic  or  mixed-linear 
parabolic  equations.  Furthermore,  the  oxygen  pressure  dependence  is 


Neither  the  oxidation  kinetics  nor  the  oxygen  pressure  dependence  are  satis¬ 
factorily  described  by  present  models  of  the  oxidation  mechanism  of  silicon. 
Alternative  models  based  on  simple  thermal  lattice  disorder  or  nonstoichio- 
metric  defects  were  formulated  and  also  found  to  be  unsatisfactory. 


Analyses  of  available  electrical  conductivity  data  for  thermally  grown 
silica  show  discordant  values  for  ion  and  electron  conduction. 


It  is  recommended  that 


1.  Oxidation  kinetics  in  dry  (<  1  ppm  H2O)  be  obtained  as  a  function 
of  T  and  p with  concurrent  studies  of  the  microstructure  and  chemical 

analyses  of  the  silica. 

2.  The  total,  ionic  and  electronic  conductivity  be  obtained  on  silica 
using  cells  of  configurations  listed  in  the  text.  Particular  attention 
should  be  directed  to  low  resistivity,  ohmic  contacts  and  to  maintaining 
the  I-V  variables  below  dielectric  breakdown  strength  of  SiO^. 


+  See  list  of  references  on  page  12. 
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